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Abstract 

We propose a method to study the degeneracy of a trapped atomic 
gas of fermions through the relaxation of the motion of a test particle. In 
the degenerate regime, and for an energy of the test particle well below 
the Fermi energy, we show that the Fermi-Dirac statistics is responsible 
for a strong decrease of the relaxation rate. This method can be used to 
directly measure the temperature of the fermionic gas. 

PACS numbers: 05.30.Fk, 03.75.Fi, 32.80.Pj 

New techniques for cooling atomic gases, based either on laser manipulation 
or evaporative cooling, have led to spectacular progresses in the realization of 
degenerate atomic samples. The most striking example is the Bose-Einstein 
condensation of alkali vapors ^ |^ and more recently, of an atomic hydrogen 
gas These techniques can also be applied to fermionic samples, in which 
spectacular phenomena such as Cooper pairing of atoms |^ , or inhibition of the 
spontaneous emission of an excited atom 1^ have been predicted. In this article, 
we present a simple and powerful tool to study the degeneracy parameters of 
such a fermionic sample, which we consider ideal for simplicity. 

Numerous possible ways to analyze an atomic Fermi gas in the degenerate 
regime have been recently studied from an theoretical point of view. Several 
authors focus on the interaction of the atomic cloud with light: modifications of 
the refraction index or of the absorption coefficient of the gas [0, ||] , reduction 
of the spontaneous emission rate of an excited atom inside the atomic cloud 
and angular dependence of the radiation pattern |T^. Spectacular effects 
are expected when the Fermi momentum pf of the cloud is larger than the 
photon momentum. Unfortunately in this regime, a quantitative analysis of the 
interaction of the gas with resonant light is difficult. Indeed it corresponds to 
nX^ > 1, where n is the spatial density of the gas and A is the wavelength of the 
light. Light propagation in the medium is then strongly affected by multiple 
scattering effects, which are difficult to handle. 

Other propose diagnostics of Fermi degeneracy of a trapped gas involve the 
study of its spatial distribution, either for a pure fermionic sample ||ll| , 

or for a mixture with a Bose-Einstein condensate These methods 

are well adapted to study the region where the temperature T is around the 
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Fermi temperature Tp, but become less sensitive to temperature in the strongly 
degenerate regime where T <^ Tp. Also, as pointed out in [|lOj, one can take 
advantage of inelastic processes inside the gas to collect information about its 
degeneracy parameter. 

In the present paper we investigate the dynamics of a probe particle (P) in 
the same trapping potential as the Fermi gas. We assume that we can prepare 
P on an arbitrary orbit of the trapping potential. An experimental procedure 
for such a preparation is outlined in the last part of this paper. Neglecting any 
inelastic process, the only way to change the trajectory of P is an interaction 
with the Fermi gas. We show that the coUisional dynamics of P gives access to 
the statistical properties of the Fermi gas, in particular its temperature T. 

The trapping potential for the test particle C/(r) is assumed to be harmonic 
and isotropic, with a frequency aj/27r. For simplicity we neglect the modification 
to the potential due to the mean field interaction between P and the Fermi 
gas. We suppose also that the same trapping potential acts on the Fermi gas. 
We calculate the dynamics of P in the Fermi sea by means of the Boltzmann 
transport equation n& in the local density approximation: 



9w „ „ „ dw 



(1) 



where 'w{r,p,t) is the phase space density of p at time t. We put v = p/M, 
where M is the mass of P. With a notation similar to |l6j the coUisional contri- 
bution to the quantum Boltzmann transport equation for P is: 



dw 
'dt 



(r,p,t) = -^ / d%d^n {wf{l~n-w'ni-f)) |v-v^| 



(2) 

where a is the coUisional cross section for interactions between P and fermions. 
Collisions occur in the low energy range where s-wave scattering is dominant 
hence the cross-section is isotropic and independent of energy. The first part of 
the coUisional integral corresponds to the decay of the phase space density in 
r, p due to a collision between P and a fermion with momentum p/: p + p/ — > 
p' + p' f- The final relative velocity is pointing in the direction given by fl. The 
second part of the integral describes the reverse process, and we put Vf=pf /m, 
where m is the fermionic mass. We assume that a single probe particle is present, 
or that the gas of probe particles is sufficiently dilute to be treated as an ideal 
Boltzmann gas. In the latter case we suppose that the number of test particles is 
sufficiently low that the relaxation does not significantly perturb the distribution 
of the Fermi gas. We use the abbreviations w' — w{r,p',t), f — /^(r, p/) and 
/' — fT{r,p' f). The quantity /T(r,p) represents the steady-state phase space 
distribution of the Fermi gas in the local density approximation: 

" 1 -I- exp ((p2/2m + Uir) - n) /(ksT)) 

where is the chemical potential. The Pauli exclusion principle is represented in 
the coUisional integral by the factors 1 — /' and 1 — / , which give the occupation 
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of the final state of the collision . The local density approximation is valid 
when tfPf ^ ?i where rp and pp are the sizes in position and momentum space 
of the Fermi gas. For a temperature much lower than the Fermi temperature 
Tp (where fc^Tp = {QN)^/'^ huj), this requires 3> 1. Indeed, one gets in this 
case rp = (487V)i/^aHO andp;^ = (48A^)^/*^pho, where oho = {h/{muj)y/'^ and 
Pao = {rnhujy^"^ are the spatial and momentum extensions of the ground state 
harmonic oscillator. 

We shall assume that the u;(r, p) is initially a distribution centered in rg, po, 
much narrower than the steady-state distribution of (||): Wss oc exp(— (p^/2M + 
U{r))/{kBT)). As the collisional relaxation proceeds the population of the nar- 
row peak is transferred to a broad distribution proportional to Wss- In the 
following we focus on the initial stage of this relaxation phenomenon, namely 
the decay of the narrow peak, which occurs at the rate Tt{vo, Po) deduced from 

rT(ro, Po) - ^ y cPn /(I - /') |vo - v;| (4) 

This rate could also be derived from the Fermi- Golden-Rule, within the local 
density approximation. 

We now consider three different classes of trajectories for P and discuss how 
the damping of those trajectories is affected by the Fermi statistics of the cloud. 
Most of the calculations assume equal mass for the P and the fcrmionic atoms. 
Such a condition is approximately realized if P and the fermions are isotopes of 
the same element {i.e. ^Li and ^Li, '^^K and '*°K). For each trajectory class we 
calculate the multiple integral for an arbitrary temperature numerically and 
we derive scaling laws for interesting limiting cases. 

The first situation consists of P at rest in the trap center. For a Boltzmann 
gas with the same number of atoms, we would expect r^°''-'(0,0) = navth, 
where n is the spatial density at the center of the trap for this gas and wth 
is the most probable speed {vth — i&kBT/ {nmY)^^'^). To put in evidence the 
effects of the Fermi statistics, we plot in FigJ^ the rate FtCO, 0) normalized 
by r^°''''(0, 0). At high temperature {T > Tp) the deviations due to Fermi 
statistics are negligible. On the other hand, for T < Tp, these deviations are 
spectacular and we find that Ft/F^"' ' cx T"^. This power law dependence 
originates from two different phenomena both related to Fermi statistics, (i) 
For the Boltzmann gas, n cx T~^^'^ and v cx VT, while in the fermionic gas 
the spatial density and the most probable speed remain constant for vanishing 
T. This gives account for a factor cx T. (ii) For T ^ Tp, P has a finite 
collision probability only with fermions within an energy interval AE ^ ksT 
(see Fig||)at the surface of the Fermi sphere, The solid angle Ail available 
for the allowed fermionic final states is also proportional to T . Therefore the 
collisional rate is reduced by an additional factor AE Af2 cx . 

This situation is well suited for determining the temperature of the Fermi 
gas in the degenerate regime. Its absolute calibration depends on the exact mass 
ratio between P and the fermions. Consequently in Fig.|l| we also plot rT(0,0) 
for the specific case of a cesium atom (^^'^Cs) as P for a lithium gas (^Li). 
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We now consider a second type of trajectory consisting of an oscillation of P 
with an energy E, and a angular momentum L which we set equal to zero (see 
insert of Figjsj). To calculate the rate at which P is ejected from this trajectory 
by collisions, we suppose that the relaxation is slow with respect to the period 
of an oscillation (collisionless regime). Since rT(r, p) is not constant over the 
trajectory, we define an average collision rate: 

jT{E,L = Q) = ^(j^TT{r{t),p{t))dt . (5) 

Depending on the energy of the oscillation and the temperature of the fermionic 
gas, we find different regimes (see Fig.||). As expected, for T > Tp the relaxation 
rate is the same as the one predicted for a Boltzmann gas. It doesn't depend 
much on the energy of the excitation even for E > ksT (as long as we assume 
a constant s-wave elastic cross section). Indeed, when E grows the relative 
velocity between P and the cloud increases as E^^^ while the fraction of the time 
P spends within the cloud decreases as leading to a constant average rate. 

For the degenerate case T < Tp, three energy domains have to be considered for 
P. (i) For E < UbT < Ep we recover the rate rT(0,0) displayed in Fig.|l|. (ii) 
For ksT < E < Ep, the rate varies as E'^. This can be easily understood at zero 
temperature with E <C Ep. In this case only collisions with fermionic particles 
with energy close to Ep contribute (first factor E), and the solid angle available 
for the allowed final states brings an extra factor E. (iii) For /cbT < Ep < E, 
all final states for the fermions after the collision lie above the Fermi surface so 
that the inhibition due to statistics in no longer effective. One recovers in this 
case a rate independent from E as for a Boltzmann gas. 

The last type of trajectory considered in this paper consists of a circular orbit 
of the test particle in the trapping potential (L — E/uj). The corresponding 
decay rate is plotted in Fig.^ as a function of the mechanical energy of P. For 
a weakly degenerate gas we recover the same result as for a Boltzmann gas. 
When one increases the energy of P, the damping rate is constant up to UbT 
and then decreases, as P is orbiting outside the cloud. In the degenerate case 
the damping rate presents a resonant behavior around E = Ep. For E <^ Ep 
this rate is decreased because of Pauli's exclusion principle, while for E ^ Ep 
it is small since P is outside the Fermi cloud. 

We now address the preparation of P with arbitrary initial position and 
momentum. The simplest idea is to exploit the difference in the ground state 
hyperfine splitting between the various alkali atoms, or between two isotopes of 
the same species. Consider for example the specific case of Lithium atoms. If 
one starts with a mixture of ^Li and ^Li, one can perform radiofrequency evapo- 
ration around a frequency t'r.f. ^ 804 MHz, corresponding to the ^Li hyperfine 
splitting. One can prepare in this way an ultracold sample of ^Li, playing the 
role of P, without eliminating any fermionic atom ^Li, whose hyperfine splitting 
corresponds to 228 MHz The state obtained in this way corresponds to the 
first situation considered in this paper. One can then prepare P on an arbitrary 
trajectory by means of successive optical stimulated Raman transitions | pO| . 
Due to the isotopic shift of the Li resonance line (10 GHz), these transitions can 
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be made isotopically selective. In a realistic case, we can consider 10® fermions 
in an isotropic magnetic trap with lo jl-K = 100 Hz (Tp — 4 /xK, rp = 0.17 mm). 
The density of the fermionic gas is 5 10^^ cni~^, giving a mean field energy 
created by the Fermi cloud on P equal to 10 nK |2^, |2^, which is negligible 
with respect to Ep, as assumed in this paper. 

To summarize we have shown in this paper that the coUisional relaxation of 
a probe particle imbedded in a Fermi atomic cloud gives a direct access to the 
quantum degeneracy of this gas. In the degenerate regime P can be regarded as 
an excitation "frozen" by Pauli's exclusion principle. One can determine both 
the temperature of the Fermionic cloud from the value of the relaxation rate F 
for E ^ 0, and the value of the Fermi energy exploiting the resonant behavior of 
F for E ~ Ep. For simplicity we have considered here a non-interacting Fermi 
gas, but it is clear that this method can be extended to study the effects of 
interactions onto the fermionic excitation spectrum. In particular we plan to 
address the consequences of Cooper pairing of the fermions ||, |3|, |2j, in a 
subsequent paper. 
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Figure 1: CoUisional rate rT(r = 0, p = 0) of a probe particle (P) of mass M at 
rest in the bottom of a trap containing a fermionic gas of temperature T (atomic 
mass m). The temperature is plotted in units of the Fermi temperature Tp, and 
Ft in units of the rate F^°'- for a Boltzmann gas at the same temperature and 
with the same number of atoms. Squares: M = m, circles: M = (133/6)m 
(case of a ^Li gas probed by ^^"^Cs). 



Figure 2: Collision between a fermion and P at rest. The fermion has an initial 

momcntiim equal to the Fermi momentum pp. p'p and p'^ arc the momenta of 
P and the fermion after collision. The circumference with diameter pp passing 
through the centre of the Fermi sphere gives the possible final states in the 
case of equal mass. The shell of thickness mkBT/pF represents the final states 
available to the Fermi particle. 



Figure 3: Damping rate Jt{E, L = 0) for a linear oscillation of P as a function 
of its excitation energy E. E is expressed in units of the Fermi energy and 7t 
is normalized by the rate navp, where n is the fermionic density at the bottom 
of the trap and vp the Fermi velocity. Circles: T — Tp, squares: T — Tp/4, 
diamonds: T = Tp/16, triangles: T = Tp/QA. The inset shows a trajectory of 
P through the fermionic cloud. 



Figure 4: Damping rate ')t{E,L = E/u) for a circular orbit of the TP as a 
function of its excitation energy E. The normalizations and symbols are the 
same as in Fig.3. 
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